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FOREWORD 

In  this  report  there  is  offered  an  improved  method  for  the  fairing  of  ship  lines,  for 
which  a need  has  been  indicated  by  experience  with  a computation  system  for  surface 
wave  trains.  Illustrative  ship  plans  for  tankers  have  been  obtained  from  Mr  Riesenberger 
of  the  Maritime  Administration  and  for  warships  from  Mr  Keane  of  the  Naval  Ship 
Engineering  Center,  The  manuscript  for  the  report  was  completed  by  15  July  1977 
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ABSTRACT 

Straight  lines,  circles,  and  splines  are  used  in  most  mathematical  representations 
of  ship  lines.  The  representations  are  only  piecewise  analytic.  Discontinuities  in 
curvature  occur  along  waterlines  and  section  lines  The  discontinuities  upset 
computations  of  velocity  distribution  in  the  flow  around  the  ship  The  discontinuities 
are  smoothed  when  the  lines  are  approximated  by  orthonormal  polynomials  The 
discontinuities  are  eliminated  when  the  lines  are  simulated  with  the  aid  of  radicals 
A smooth  simulation  of  the  Series  60.  Block  0 60  ship  model  is  given  by  a radical 
representation 


INTRODUCTION 


The  traditional  representation  of  a body  of  revolution  is  a table  of  radial  distances 
from  the  axis  of  symmetry  at  a series  of  stations  along  the  axis  of  symmetry  The 
traditional  representation  of  a foil  or  a ship  is  a table  of  offsets  at  a rectangular  grid 
on  the  median  plane  of  the  foil  or  the  ship  The  offsets  for  the  ship  define  the  waterlines 
and  the  section  lines  on  the  surface  of  the  ship 

An  introduction  to  the  mathematical  representation  of  the  surface  of  a ship  has 
been  published  in  a book  by  Kuo1  Included  in  the  book  are  analyses  of  the  curved 
sections  of  the  ship  from  stem  to  stern  Reviews  are  given  for  some  published  methods 
for  the  fairing  of  ship  lines  Further  analyses  of  ship  lines  have  appeared  in  more 
recent  publications2  9 

The  flow  over  a body  has  been  simulated  by  Hess  and  Smith12  14  with  a source 
distribution  over  a polyhedron  whose  vertices  he  on  the  surface  of  the  body  The 
normal  and  the  area  at  the  surface  are  assumed  to  be  the  normal  and  the  area  at  a 
face  of  the  polyhedron  The  source  density  is  assumed  to  be  uniform  over  each  face 
of  the  polyhedron,  and  the  velocity  in  an  ideal  flow  over  the  polyhedron  is  taken  to 
be  the  velocity  at  the  midpoint  of  each  face  Accurate  agreement  between  the-  polyhedral 
simulation  and  an  analytical  representation  is  possible  for  a polyhedron  with  a large 
enough  number  of  faces  The  velocity  for  potential  flow  over  a real  polyhedron  would 
be  infinite  at  each  edge  of  the  polyhedron  The  flow  over  a dihedral  angle  is  analysed 
in  Appendix  C 

For  an  irregular  shape  it  would  be  necessary  to  use  .1  polyhedral  simulation,  but 
for  a smooth  shape  it  should  be  possible  to  take  advantage  of  the  smoothness  in  order 
to  increase  accuracy  and  reduce  cost  The  interaction  between  two  points  in  an 
unbounded  fluid  is  given  by  the  simple  inverse  square  law.  whereas  the  interaction 
under  a free  surface  is  given  by  the  derivatives  of  an  expensive  Havelock  integral  A 
reduction  of  cost  by  the  use  of  fewer  data  is  vital  for  surface  ships 

It  has  been  the  goal  of  naval  architects  to  achieve  smooth  waterhnes  Smoothness 
was  achieved  by  eye  The  lines  of  a ship  were  laid  out  in  a mold  loft,  where  ducks* 
and  splines  were  used  in  fairing  the  lines  A duck  is  a lead  weight  with  a protruding 
tongue,  and  a spline  is  a slender  strip  with  a groove  along  the  edge  The  tongue  of 
the  duck  engages  the  groove  of  the  spline  and  holds  it  in  place  The  tongue  of  the 
duck  is  finite  in  breadth  Translation  and  rotation  of  the  duck  controls  the  position 
and  the  slope  of  the  spline  The  curvature  is  not  continuous  across  the  node  where 
the  duck  is  in  contact  with  the  spline 

Much  attention  has  been  devoted  to  mathematical  splines  in  the  literature  An 
analysis  is  to  be  found  in  Hildebrand's  Introduction  to  Xumcrica/  Analysis"  The  length 
of  a curve  is  subdivided  into  intervals  The  ordinate  in  each  interval  is  expressed  in 
terms  of  the  abscissa  by  a cubic  polynomial  Inasmuch  as  the  cubic  polynomial  has 
four  coefficients,  it  can  be  made  to  meet  four  conditions  It  can  be  made  to  have 
specified  ordinates  at  the  two  ends  of  the  interval  It  can  be  made  to  maintain 
continuity  of  first  and  second  derivatives  across  boundaries  between  intervals  The 
coefficients  of  the  spline  approximations  are  the  solution  of  a tridiagonal  system  of 
equations  Recurrence  expresses  each  successive  approximation  in  terms  of  a pair  of 
parameters  which  are  determined  by  conditions  at  the  ends  of  the  curve 


Vspime  weights  are  called  dolphins  in  Webster's  unabridged  dictionary' 
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Linearized  cubic  splines  have  been  analysed  by  Mehlum8  From  the  calculus  of 
variations  he  has  shown  that  the  minimization  of  the  square  of  the  second  derivative 
leads  to  the  cubic  spline  The  elastic  curve  is  the  equilibrium  configuration  of  a physical 
spline.  The  minimization  of  the  square  of  the  curvature  along  an  elastic  curve  leads 
to  the  special  case  of  a physical  spline  to  which  forces  are  applied  but  not  torques 
An  analysis  of  the  physical  spline  is  given  in  Appendix  D 

Minimization  of  the  square  of  the  curvature  has  been  used  for  smoothing  in  the 
AUTOKON  system  by  Mehlum  For  the  representation  of  smooth  lines  he  has  used 
sequences  of  circular  arcs.  Such  representations  have  application  in  the  fabrication 
of  ships  where  flame  cutters  are  guided  along  straight  lines  or  along  circles  Let  the 
successive  points  along  a curve  be  joined  by  chords  which  span  the  intervals  betiveen 
the  points  Continuity  of  ordinates  from  interval  to  interval  is  achieved  if  the  chords 
are  spanned  by  circles  with  centers  on  the  perpendicular  bisectors  of  the  chords 
The  circular  arc  between  two  points  and  the  radii  to  the  ends  of  the  arc  are  the 
boundaries  of  a sector  of  a circle.  Continuity  of  slopes  from  interval  to  interval  is 
achieved  if  the  centers  of  the  circles  are  located  at  the  intersections  between  boundary 
radii  and  perpendicular  bisectors  The  representation  of  a curve  as  a sequence  of 
circular  arcs  can  be  constructed  in  a straightforward  progression  in  terms  of  the 
radius  of  the  first  circle  The  construction  of  a smooth  curve  requires  an  adjustment 
of  the  first  radius  by  iteration  That  the  AUTOKON  system  has  been  successful  for  the 
fairing  of  150  ships  has  been  reported  by  Mehlum  That  special  skill  by  an  operator 
is  required  in  the  fairing  process  has  been  reported  by  Doornbos9 

Any  simulation  with  the  aid  of  splines  can  be  only  piecewise  analytic  Furthermore, 
freedom  from  the  limitations  of  traditional  methods  is  possible  now  that  plotters  and 
machine  tools  can  be  driven  by  servo  motors  under  the  control  of  computers 

A conformal  mapping  of  a section  line  on  a circle  has  been  proposed  by  von  Kerczek 
and  Tuck6  Their  representation  of  a section  line  is  just  an  expression  of  Cartesian 
coordinates  in  terms  of  a parametric  angle  by  trigonometric  polynomials  Inasmuch 
as  the  representation  is  horizontal  at  the  base  line  and  is  vertical  at  the  waterline, 
it  is  not  suitable  for  most  section  lines. 

A mathematical  definition  of  a surface  is  obtained  when  Cartesian  coordinates  of 
the  surface  are  expressed  as  parametric  functions  of  surface  coordinates  The 
representation  of  Cartesian  coordinates  in  terms  of  surface  coordinates  makes  possible 
the  computation  of  normals  and  areas.  A quadrilateral  boundary  on  the  curved  surface 
can  be  mapped  into  a square  boundary  on  a plane  map  The  connection  between 
Cartesian  coordinates  on  the  quadrilateral  and  the  surface  coordinates  m the  map  is 
indeterminate  A point  on  the  surface  remains  on  the  surface  for  any  arbitrary 
transformation  which  moves  the  point  parallel  to  the  surface 

A study  of  orthogonal  transformations  has  been  made  by  Miloh  and  Patel7  They 
proposed  trajectory  integration  to  trace  lines  which  are  orthogonal  to  section  lines 
They  considered  also  an  identification  of  coordinate  lines  with  traces  of  minimum 
and  maximum  curvature  Such  an  identification  would  give  a unique  mapping 
transformation,  but  they  found  that  the  coordinate  lines  would  have  awkward 
curvatures 

The  connection  between  Cartesian  coordinates  and  map  coordinates  is  determinate 
when  the  map  transformation  is  isometric  The  Cartesian  coordinates  are  expressed 
in  terms  of  distance  along  lines  in  a network  on  the  surface  An  isometric  mapping 
was  used  already  in  1960  for  the  Taylor  Standard  Ship17  The  use  of  distance  along  a 
curve  as  a parametric  variable  has  been  proposed  again  recently  by  Jancaitis  and 
Junkins16 
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The  parametric  representations  can  be  expressed  with  the  aid  of  polynomial 
approximations.  Past  experience  with  higher  order  polynomials  has  shown  that  the\ 
are  difficult  to  control  when  they  are  based  upon  data  which  have  equally  spaced 
arguments.  They  are  not  difficult  to  control  when  they  are  based  upon  data  which 
have  arguments  with  the  spacing  of  the  roots  of  a Chebyshev  polynomial  Even  when 
the  arguments  are  equally  spaced,  the  polynomials  are  easy  to  control  when  least 
squares  are  applied  to  a large  number  of  closely  spaced  data 

A mathematical  representation  of  the  Series  60.  Block  0 60  ship  model  was  reported 
at  the  First  International  Conference  on  Numerical  Ship  Hydrodynamics18  Data  on  the 
hull  configuration  were  derived  from  a table  of  offsets  in  Todd’s  report15  The  hull  has 
flat  sides  and  a flat  bottom  which  are  joined  by  a curved  surface  The  flat  bottom 
intersects  the  curved  surface  along  a line  of  tangency  The  actual  models  had  various 
stern  configurations  with  propellers  and  rudders  Only  the  simplest  configurations  at 
the  stem  and  at  the  stern  are  within  the  scope  of  the  present  project  The  actual 
models  had  a length  of  20  ft  from  fore  perpendicular  to  aft  perpendicular  The 
mathematical  model  has  been  scaled  to  have  the  same  dimensions 

The  tabular  data  in  Todd's  report  are  too  sparse  to  define  the  configuration 
mathematically  An  attempt  to  force  the  section  lines  to  be  tangent  to  the  flat  bottom 
on  the  tabulated  line  of  tangency  caused  irreducible  wiggles  in  a few  of  the  section 
lines  It  was  necessary  to  assume  that  each  section  line  intersects  the  bottom  at  a 
finite  angle  The  table  of  offsets  was  augmented  with  additional  section  lines  and 
additional  waterlines  near  the  edges  of  the  sides  of  the  model 

Each  side  of  the  hull  was  simulated  by  a single  quadrilateral  which  extends  from 
the  fore  perpendicular  to  the  aft  perpendicular.  The  arguments  of  polynomial 
approximations  w'ere  proportional  to  the  longitudinal  distance  to  a cross  section,  and 
to  the  lateral  distance  along  the  section  line.  The  table  of  offsets  was  adjusted  by  trial 
until  the  waterlines  and  the  section  lines  were  fair  in  CalComp  plots 

Actual  waterlines  have  discontinuities  in  curvature  where  the  forebody  joins  the 
midbody  and  where  the  midbody  joins  the  aftbody  Actual  section  lines  have 
discontinuities  in  curvature  at  the  fore  perpendicular  where  the  straight  stem  bends 
outward,  at  the  midship  section  where  the  curved  bilge  joins  the  straight  side,  and 
at  the  aft  perpendicular  where  the  straight  stern  bulges  outward  Section  lines  near 
the  stern  are  flat  at  the  bottom,  convex  at  the  bilge,  concave  at  the  bend,  and  convex 
at  the  bulge  Associated  with  discontinuities  in  curvature  are  discontinuities  in  velocity 
which  would  disrupt  a numerical  analysis  of  velocity  Furthermore,  the  mathematical 
model  had  a dihedral  angle  where  the  sides  met  the  bottom  Any  potential  flow  over 
a convex  dihedral  angle  would  have  infinite  velocity  at  the  vertex  of  the  dihedral  angle 
This  might  explain  a bizarre  velocity  distribution  which  was  obtained  near  the  junction 
between  the  sides  and  the  bottom 

A more  realistic  flow  field  would  have  separation  with  the  origination  of  vortex 
sheets  from  the  vertex  of  each  dihedral  angle  The  analysis  of  flow  with  vortex  sheets 
is  beyond  the  scope  of  the  present  project. 

There  is  clearly  a need  for  a truly  mathematical  formulation  which  provides  all  of 
the  characteristic  features  of  a flat  bottomed  merchant  ship  Section  lines  must  be 
tangent  to  the  flat  bottom  There  must  be  provision  for  a nearly  parallel  midbody 
There  must  be  the  right  number  of  inflection  points  and  regions  of  maximum  curvature 
in  both  waterlines  and  section  lines  It  is  the  objective  of  this  report  to  offer  a 
mathematical  formulation  which  is  suitable  for  global  flow  over  a ship's  hull  Not 
included  in  the  report  are  formulations  for  local  flows  over  appendages  such  as  domes, 
propellers,  or  rudders 
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There  are  two  ways  to  functionalize  a ship  hull  In  one  way  the  section  lines  are 
expressed  by  simple  formulae  with  parameters  which  vary  lengthwise  In  the  other 
way  the  waterlines  are  expressed  by  simple  formulae  with  parameters  which  vary  with 
elevation.  The  second  way  has  been  found  to  be  more  satisfactory,  provided  seven 
equally  spaced  sections  were  expressed  with  simple  functions.  This  was  enough  to 
establish  waterlines  at  all  elevations  and  section  lines  at  all  stations. 

Two  types  of  parameters  can  be  recognized  One  type  is  only  weakly  determinate, 
and  can  be  expressed  therefore  by  simple  formulations.  The  other  type  is  strongly 
determinate,  and  must  be  expressed  as  the  solution  of  simultaneous  equations 

In  the  present  report  there  are  documentations  of  programs  and  subroutines  which 
are  specific  to  ship  lines  In  other  reports  there  will  be  documentations19  20  of  subroutines 
which  are  of  general  use  for  polynomials  and  matrices 

CARTESIAN  COORDINATES 

In  the  simulation  of  offsets  at  the  surface  of  the  ship  it  is  convenient  to  use  a 
special  coordinate  system  in  which  x is  the  distance  backward  from  the  fore 
perpendicular,  y is  the  offset  from  the  median  plane,  and  z is  the  distance  upward 
from  the  flat  bottom  Associated  with  the  Cartesian  coordinates  are  unit  vectors  i.j.k 
in  the  directions  of  increasing  coordinates  x.  y.  z 

WATERLINES 

In  accordance  with  a fundamental  philosophy  of  ship  design,  waterlines  and  section 
lines  are  smooth  curves  with  asymptotes  which  are  parallel  to  the  median  plane  of 
the  ship  Smoothness  may  be  achieved  with  simple  mathematical  formulations.  The 
simplest  curve  with  straight  asymptotes  is  the  hyperbola,  in  which  case  the  ordinates 
are  expressed  in  terms  of  the  abscissae  by  radicals  The  ordinate  y for  a simple 
hyperbola  is  expressed  in  terms  of  the  abscissa  x by  the  equation 

y = a \ pz  + az  ( 1 ) 

where  the  parameter  n is  defined  by  the  equation 

<T  = x - x0  (2) 

while  x0  is  the  locus  of  the  center  of  the  hyperbola,  of?  is  the  semidiameter,  and  a 
is  the  slope  of  the  asymptotes 

It  is  possible  to  string  together  a set  of  hyperbolas  with  different  positions  and 
dimensions  to  form  excellent  simulations  of  conventional  waterlines  Let  a waterline 
be  represented  by  the  equation 

y a0  * o,\  /S,2  + it,2  + a2\  P2Z  + o2  + a3\  P3  - o32  • o4\  pt 2 • a42  (3) 

where  the  parameters  are  defined  by  the  equations 
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The  constant  n0  controls  the  amount  by  which  the  waterline  is  displaced  off  the  median 
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plane,  and  the  coefficients  at , a2.  a3.  a4  control  the  amplitudes  by  which  the  waterline 
swings  in  and  out.  The  asymptotes  of  the  curves  which  are  formed  by  pairs  of  hyperbolas 
are  parallel  to  the  median  plane  when  the  coefficients  are  related  in  accordance  with 
the  equations 

a2  = - a,  a4  = — n3  (8) 

The  parameters  O0.aj.a3  can  be  computed  by  the  solution  of  three  equations  which 
constrain  the  offsets  to  have  specified  values  at  the  fore  perpendicular,  at  the  midship 
section,  and  at  the  aft  perpendicular 

The  constants  pt,  pz.  p3,  /?4  control  the  curvature  of  the  waterline  in  the  four  regions 
where  the  curvature  is  a maximum.  The  curvature  has  fore  and  aft  symmetry  when 
the  constants  are  related  in  accordance  with  the  equations 

Pi  = 04  0Z  = 03  O) 

The  constants  have  so  little  effect  on  the  waterline  that  they  are  almost  arbitrary 
They  are  adjusted  by  trial  to  give  good  simulations  of  conventional  waterlines. 

The  constants  xl,x2,x3.xi  determine  the  position  and  the  slope  of  the  waterline 
near  two  points  of  inflection.  The  constants  have  a strong  effect  on  the  waterline  They 
can  be  adjusted  by  trial  until  the  offsets  are  equal  to  specified  values  at  four  specified 
positions  along  the  length  of  the  ship. 

An  attempt  was  made  to  use  least  squares  for  the  adjustment  of  constants  to  a 
seven- station  data  set.  There  was  convergence  when  the  data  were  equally  spaced, 
but  there  was  divergence  when  the  data  were  moved  too  close  to  the  fore  perpendicular 
and  the  aft  perpendicular  The  rate  of  convergence  or  divergence  was  increased  when 
sequential  iteration  of  the  four  constants  was  replaced  by  the  simultaneous  evaluation 
of  the  constants  through  inversion  of  a 4 * 4 matrix.  Thus  there  are  limitations  on 
the  range  of  application  of  the  radical  representation 

In  a series  of  experiments  the  waterlines  were  simulated  with  the  product  of  sums 
of  radicals,  but  this  complication  was  not  superior  to  the  simple  sum  of  four  radicals 


SECTION  LINES 


The  bend  in  a section  line  is  expressed  by  a linear  combination  of  a straight  line 
and  a hyperbola.  It  was  found  possible  to  orient  line  and  hyperbola  so  that  one  leg 
of  the  combination  was  parallel  to  the  median  plane  The  bulge  in  a section  line  is 
achieved  through  t.h<  division  by  another  radical.  If  the  constants  in  both  radicals  are 
the  same,  then  the  curve  is  expressed  by  the  sum  of  a constant  and  the  derivative 
of  a radical  as  illustrated  by  the  equation 


Pa  + \ 1 1-  p~a‘- 
\ 1 • pzaz 


\ 1 + P'-a" 


where  the  parameter  a is  defined  by  the  equation 

a z a 


( 11) 


and  a,  p are  constants  The  curve  which  is  expressed  by  these  functions  is  S - shaped 
with  horizontal  asymptotes  Termination  of  the  curve  at  the  bottom  where  2 0 is 

achieved  if  the  functions  are  multiplied  by  the  factor 


n 


where  e is  a constant  This  factor  vanishes  like  the  square  root  of  2 as  z - 0,  and 
approaches  unity  as  z ->  The  curve  which  is  expressed  by  the  factor  is  U— shaped 
with  parallel  asymptotes  The  bottom  of  the  curve  osculates  a circle  if  the  radius  of 
curvature  is  equal  to  unity.  Substitution  of  the  factor  into  Equation  (6)  in  Appendix  D 
shows  that  the  radius  of  curvature  is  unity  if  e = | An  opening  at  the  end  of  waterhnes 
or  section  lines  can  be  closed  with  the  simulation  of  a circle  if  e is  scaled  to  match 
the  width  of  the  opening 

It  is  possible  to  combine  functions  to  form  excellent  simulations  of  conventional 
section  lines.  Let  a section  line  be  represented  by  the  equation 

where  the  parameter  a is  defined  by  the  equation 


The  constant  a determines  the  position  of  the  bend  The  constants  fi.  y determine  the 
angle  between  asymptotes  of  the  bend  The  constant  <5  dete  .lines  the  curvature  of 
the  bulge.  The  constant  € determines  the  curvature  of  the  bilge  The  constants  .4,  B.  C 
can  be  computed  by  the  solution  of  three  equations  which  constrain  the  offsets  to 
have  specified  values  at  three  elevations. 


DERIVATIVES 

The  computation  of  areas  and  normals  requires  the  evaluation  of  partial  derivatives. 
The  offset  y at  the  surface  is  expressed  as  a function  of  coordinates  x,  z on  the  median 
plane  by  the  equation 

y = f(x,  z)  (15) 

The  function  is  given  by  Equation  (3)  where  the  parameters  ct,,ct2.ct3,ct+  are  the  only 
functions  of  x,  and  the  other  parameters  are  functions  of  z Differentiation  with  respect 
to  x is  expressed  by  the  equation 

9/  a | it,  n ,a  , a3cr3  | a4a4 

dx  + a,2  + azz  V/g32  + a3z  v'/?42  + a} 

while  differentiation  with  respect  to  z is  expressed  by  the  equation 

df  da0 
dz  dz 


v'/9,2  + a 


of  sever 


2 d“l 

+ V /Sa2 

+ a/ 

da  2 

v P3  + a 

, dn  , 

3 i 

• 

0. 

1 dz 

dz 

dz 

iZ 

dfit 

<*2 

Pz 

dp 2 4 

a 3 Ps 

dP} 

a4P4 

dP4 

Z dz 

+ a/ 

dz 

\ + a 

32  dz 

" N P4^ 

dz 

rfx, 

f*2 

°z 

dx2 

0.3  <7.1 

dx3 

Otg4 

dx4 

,2 

+ az 

dz 

\ P3Z  + a 

? dz 

V P/  ■ 

dz 

• parameters 

can 

be  det 
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Differentiation  in  Equation 


(13)  leads  to 


the  equation 


dy 

dz 


\P 


( y 2 - 


( 1 + 6zoz)z  L 

(So  + v/l 


v 1 + yzaz 


B + C 


yzaz 


v 1 + 6zaz 


(18) 


which  can  be  used  in  the  determination  of  the  derivatives  of  the  offsets  at  the  seven 
stations 

The  vector  area  on  the  surface  per  unit  area  on  the  median  plane  is  given  by  a 
vector  cross  product  with  the  components 


df 

dx 


- — . + 1. 


V 

dz 


(19) 


Although  these  equations  are  complicated  they  must  be  used  for  the  computation  of 
area  unless  the  radical  representation  can  be  replaced  accurately  by  an  orthonormal 
polynomial  representation 


ORTHONORMAL  POLYNOMIALS 

Let  ut,  , un  be  a discrete  set  of  values  of  a parameter  u.  and  let  <f0(u).  ^m(u) 

be  a set  of  polynomials  of  progressively  increasing  degree.  The  polynomials  are 
orthonormal  if  they  satisfy  the  equations 

E VfMvM  = (20) 

i = 1 

where  6,k  is  zero  if  j k and  is  unity  if  j = k The  polynomials  are  generated  by  the 
three  term  recurrence  equation 

<Pk(u)  = ' iM  - Pk-ifk  i(u)  ~ ak-z!Pk  2<u);  (2D 

at-i ( ' 

which  is  started  with  the  polynomial 

?o(u)  = 4=  (22) 

v n 

and  uses  the  recurrence  constants  at.,  and  pk  , The  derivatives  of  the  polynomials 
are  generated  by  the  three  term  recurrence  equation 

<p'k(u)  1 -yk  t(u)  + u<p'k  ,(w)  Pk  tfk  i(w)  ak  S!fik  2(u);  (23) 

at  , ( ) 

which  is  started  with  the  derivative 

*i(u)  = 0 (24) 

The  recurrence  equations  can  be  used  to  expand  the  orthonormal  polynomials  in 
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— ^ 


ascending  powers  of  their  arguments,  but  this  expansion  is  inadvisable  because  there 
are  large  terms  of  intermediate  degree  which  cause  loss  of  accuracy  by  round  off 
Let  c0,  . cm  be  the  coefficients  in  a series  expansion  in  terms  of  the  polynomials 

An  arbitrary  function /(u)  deviates  from  ttie  series  expansion  with  a standard  deviation 
a which  is  expressed  by  the  equation 

v2  1 -J  ;/(“.)  1 <25> 

™ ,.i  V * o > 

Differentiation  with  respect  to  the  coefficients  shows  that  for  least  squares  deviation 
the  function  is  represented  by  the  equation 

n m 

/(u)  1'  1’  /(u,)^t(u,)^t(-u)  (26) 

i 1 * 0 

and  its  derivative  is  given  by  the  equation 

(i  n m 

, f(u)  /(uJ^ttuJyJtfu)  (27) 

du  t 1 t 0 

The  series  expansion  converges  to  /(it ) with  only  weak  requirements  on  the  continuity 
of  /(u) 


ISOMKTRIC  PARAN1 KTKRS 

The  mapping  of  a quadrilateral  on  the  median  plane  is  singular  at  the  bottom  where 
the  surface  is  perpendicular  to  the  median  plane  There  is  no  such  singularity  in  an 
isometric  mapping  The  Cartesian  coordinates  i.y.z  on  the  surface  are  expressed  as 
functions  of  surface  coordinates  u.  i in  t he  map  The  surface  coordinates  arc  normalized 
so  that  they  satisfy  the  inequalities 

- 1 S u S + 1 1 • v • 1 (28) 


In  the  limiting  case  of  a slender  ship  the  coordinate  x may  be  expressed  as  a linear 
function  of  the  coordinate  u.  and  the  coordinates  y.  z may  be  expressed  with 
orthonormal  polynomials  of  the  coordinate  i Discrete  values  of  r for  the  expression 
of  coordinates  are  estimated  initially  from  the  perimeter  of  a polygon  which  is  inscribed 
in  the  section  line  The  metric  l along  the  section  line  is  defined  by  the  equation 


c// 

dv 


(dz 

. \ dv  1 

\dv 

(29) 


Values  of  the  metric  are  computed  at  the  discrete  values  of  i . the  metric  is  expressed 
as  a power  polynomial  in  v,  and  is  integrated  with  respect  to  i The  discrete  values 
are  refined  by  iteration  until  the  metric  is  uniform  along  the  section  line 

A position  vector  r on  the  surface  of  the  quadrilateral  is  given  by  the  equation 


r xi  +■  yj  + ;k 


no) 


A vector  surface  element  on  the  quadrilateral  is  given  by  the  cross  product  between 
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differential  displacements  of  r for  differential  changes  in  u,  i The  vector  area  s on 
the  quadrilateral  per  unit  area  in  the  map  is  given  by  the  equation 

i j k 

s dV  d=  (3!) 

du  du  du 

dx  dy  dz 
dv  dv  dv 

The  vector  s may  be  expressed  in  terms  of  its  components  l,  m,  n by  the  equation 

s li  + ?nj  + n k (32) 

The  components  l.m.n  are  the  cofactors  in  the  determinant  of  the  Equation  (31) 
Insofar  as  an  arbitrary  function  can  be  expressed  implicitly  in  powers  of  the  surface 
coordinates,  the  function  can  be  integrated  over  the  map  Seperate  integration  with 
respect  to  each  surface  coordinate  is  made  possible  when  the  map  is  subdivided  into 
a grid  of  lines  parallel  to  the  sides  of  the  map  and  with  a spacing  proportional  to  the 
roots  of  a Chebyshev  polynomial  Then  the  matrix  of  integration  multipliers  for  the 
square  map  is  the  matrix  product  of  the  arrays  of  multipliers  for  line  integration 
The  results  of  the  smooth  simulation  are  three  matrices  A',  )',  Z which  list  surface 
coordinates  and  three  matrices  L,  \l,  .V  which  list  surface  components  for  the  grid 
points  in  the  map. 


SIMULATION 

The  waterlines  of  the  Series  60,  Block  0.60  ship  configuration  have  been  simulated 
with  the  radical  representation  of  Equation  (3),  The  parameters  iri  the  representation 
were  adjusted  by  trial  until  deviations  in  offset  were  optimum  m accordance  with  the 
Chebyshev  criterion  The  radicals  for  a representative  waterline  are  illustrated  for 
comparison  in  Figure  1 

Offsets  were  computed  at  seven  stations  along  the  length  of  the  ship  from  the 
optimum  radical  representation  The  offsets  at  the  seven  stations  were  simulated  with 
the  section  line  representation  of  Equation  (13).  The  parameters  in  the  representation 
were  adjusted  by  trial  and  by  least  squares.  Further  adjustments  were  necessary  in 
the  aft  body  in  order  to  obtain  fair  section  lines  Efficient  trial  computations  were 
made  possible  by  an  HP67  programmable  pocket  computer  The  final  parameters  in 
the  smooth  simulation  are  tabulated  in  Appendix  A 

The  deviations  in  offset  between  the  tabular  data  and  the  final  simulation  are 
illustrated  in  Figure  2.  There  is  good  agreement  near  the  bow-,  but  there  is  disagreement 
near  the  stern,  because  the  stern  post  has  been  moved  aft  to  the  aft  perpendicular 
Any  residual  discrepancy  can  be  considered  as  a local  disturbance  in  the  global 
configuration 

The  final  formulation  was  used  in  th  construction  of  a closely  spaced  table  of 
offsets,  which  is  illustrated  by  points  < Figures  3 and  4 The  table  of  offsets  was 
converted  into  an  orthonormal  polynomial  representation,  which  is  illustrated  by 
curves  in  Figures  3 and  4 There  is  accurate  agreement  between  points  and  curves 
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PROGRAMMING 


Transformations  of  ship  lines  between  tabular  formats  and  functional  formats  is 
accomplished  through  the  concatenation  of  subroutines  and  programs 

The  parameters  in  the  radical  representation  are  normalized  to  give  specific  offsets 
at  three  stations  by  the  following  subroutine. 

SUBROUTINE  NWLN  (XN,  AX,  AA,  AB,  YN) 

*****++**********+****+++*++***++***++++*»+*+******+++**************************m***m*******m 
FORTRAN  SUBROUTINE  FOR  NORMALIZATION  OF  WATERLINE 

^*^^*^m*********************************^*********** ***************************************** 

The  longitudinal  coordinates  of  three  stations  are  given  in  the  3 array  ■ The 
parameters  x1,x2,X3,x4  are  given  in  the  4 array  AX.  The  parameters  (3lrf32.p 3.^4  are 
given  in  the  4 array  AB.  The  offsets  at  the  three  stations  are  given  in  the  3- array 
YN.  New  coefficients  are  computed  by  the  solution  of  three  linear  equations  The 
coefficients  a0.  a, . a2,  a3,  a4  are  stored  in  the  5-  array  AA. 

The  offset  on  a waterline  is  computed  by  reference  to  the  following  subroutine 


SUBROUTINE  WLNO  (X,  AX,  AA,  AB,  Y) 

FORTRAN  SUBROUTINE  FOR  WATERLINE  OFFSET 

ft*.*****************.********************************* 


The  longitudinal  coordinate  x is  given  in  the  argument  X.  The  parameters  x, , x2,  x3.  x4 
are  given  in  the  4 array  AX.  The  coefficients  a0,  a,.  a2,  a3.  a4  are  given  in  the  5-  array 
AA.  The  parameters  /?, , /?a. /?3.  ff4  are  given  in  the  4 array  AB  The  offset  is  computed 
from  the  radical  representation  The  offset  y is  stored  in  the  function  - 

The  derivative  of  the  offset  on  a waterline  is  computed  by  reference  to  the  following 
subroutine. 


SUBROUTINE  WLND  (X,  AX,  AA,  AB,  D) 


FORTRAN  SUBROUTINE  FOR  DERIVATIVE  OF  WATERLINE  OFFSET 


The  longitudinal  coordinate  x is  given  in  the  argument  X.  The  parameters  x1.x2.x3,x4 
are  given  in  the  4 array  AX.  The  coefficients  a0.  a, . a2,  cx3.  a4  are  given  in  the  5-array 
AA.  The  parameters  /3, . fJ3.  are  given  in  the  4 array  AB.  The  derivative  is  computed 
from  the  radical  representation  The  derivative  dy  dx  is  stored  in  the  function 

The  parameters  in  the  section  line  representation  are  adjusted  to  give  a specific 
offset  at  one  waterline  by  the  following  subroutine 


SUBROUTINE  NSLN  (ZN,  CC,  YN) 

*********************.*********************************************% 
FORTRAN  SUBROUTINE  FOR  NORMALIZATION  OF  SECTION  LINE 
*****#*******+****************************************************% 


The  elevation  coordinate  of  the  waterline  is  given  in  the  argument  . N The  parameters 
a,(i.y.6,t.A,B.C  are  given  in  the  8 array  CC  The  offset  at  the  waterline  is  given  in 
the  argument  YN  The  parameter  C is  adjusted  to  make  the  computed  offset  equal  to 
the  given  offset 
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The  offset,  on  a section  line  is  computed  by  reference  to  the  following  subroutine 

SUBROUTINE  SLNO  (Z,  CC,  Y) 


FORTRAN  SUBROUTINE  FOR  SECTION  LINE  OFFSET 


The  elevation  coordinate  2 is  given  in  the  argument  Z.  The  parameters  a,  fi.  7,  6 , c , A . B , C 
are  given  in  the  ft  array  CC  The  offset  y is  stored  in  function  Y. 

The  derivative  of  the  offset  on  a section  line  is  computed  by  reference  to  the 
following  subroutine 

SUBROUTINE  SIND  (Z,  CC,  D) 


The  x -coordinates  of  section  lines  and  the  c coordinates  of  waterlines  are  given  in 
DATA  statements.  The  offset  data  are  given  on  TAPE^  Noise  in  the  offset  data  is  reduced 
to  least  squares  by  orthonormal  polynomial  representation  The  waterlines  are  forced 
to  pass  through  the  end  points  of  the  data  by  a linear  adjustment  Plot  instructions 
for  the  CalComp  plotters  are  recorded  on  ’APE'. 

The  section  lines  of  the  ship  are  compared  with  the  data  from  which  thev  are 
derived  by  reference  to  the  following  program 

PROGRAM  BODY  (INPUT,  OUTPUT,  TAPE1,  ’APE7) 

**#******************************************♦**********•*•••*••*•*•»•••••*•••••*•••••••••••• 

FORTRAN  PROGRAM  TO  PLOT  OFFSET  DATA  AND  SECTION  LINES 


The  x-coordinates  of  section  lines  and  the  2 coordinates  of  waterlines  are  given  in 
DATA  statements.  The  offset  data  are  given  on  TAp‘  Noise  in  the  offset  data  is  reduced 


to  least  squares  by  orthonormal  polynomial  representation  The  section  lines  are  forced 
to  pass  through  the  end  points  of  the  data  by  a linear  adjustment  Plot  instructions 
for  the  CalComp  plotters  are  recorded  on  TAPtl. 


DISCUSSION 

The  least  squares  convergence  of  a polynomial  to  a radical  is  rapid  only  within  the 
range  of  convergence  of  a Taylor  series  expansion  The  slow  convergence  of  the 
polynomial  representation  limits  its  range  of  application  The  difference  between 
representations  is  exemplified  by  their  behavior  with  increasing  arguments  The 
polynomials  increase  indefinitely  with  an  increasing  rate,  whereas  the  radicals  approach 
asymptotes  of  finite  slope 

Once  one  smooth  simulation  has  been  achieved  it  is  possible  to  generate  an  infinity 
of  other  smooth  simulations  with  the  scaling  of  breadth  and  draft  Changes  in  block 
coefficient  can  be  achieved  with  shifts  in  the  centers  of  the  hyperbolas. 

The  use  of  radicals  has  been  successful  for  the  fairing  of  a merchant  ship  with  a 
fiat  bottom  The  use  of  radicals  looks  promising  for  the  fairing  of  a complete  range 
of  ships  from  tankers  with  blunt  bows  to  warships  with  sharp  bows  Representative 
ship  plans  are  in  hand,  and  some  trial  computations  will  be  undertaken  if  time  permits 
There  will  not  be  time  in  the  present  project  for  a thorough  investigation 


CONCLUSION 

A radical  representation  has  been  successful  for  the  smooth  simulation  of  a merchant 
ship.  The  possibility  of  using  radicals  in  the  fairing  of  tankers  and  warships  should 
be  investigated  further. 
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TABLE  I 


WATERLINES  AT  SEVEN  ELEVATIONS 


WL 

*1 

*3 

*4 

0 00 

3.900000 

8 800000 

1 1 300000 

17  800000 

0.25 

0 916690 

7 708895 

12.794486 

19  0241  16 

0 50 

0 910125 

6.853267 

14  222032 

18  746849 

0.75 

0.982791 

6 375 1 24 

15.383269 

19  052597 

1.00 

0 854743 

6 049528 

16  336156 

19  942718 

1.25 

0.297610 

5.809324 

17  001219 

21  032526 

1.50 

2 000000 

5.600000 

18  000000 

22.200000 

WL 

“o 

“i 

«? 

«3 

n. 

0 00 

0.004390 

0 074161 

0 074161 

0 056375 

0 056375 

0 25 

0 142102 

0 084603 

0 084603 

0 092790 

0 092790 

0 50 

0 128860 

0 097260 

0 097260 

0 1284  10 

0 128410 

0.75 

0 148512 

0 106981 

0 106981 

0 16344 1 

0 16344 1 

1 00 

0.173393 

Oil  1567 

0 1 1 1567 

0 171716 

0 171716 

1.25 

0 144701 

0 107520 

0 107520 

0 142895 

0 142895 

1 50 

0 302918 

0.093630 

0.093630 

0 146412 

0 146412 

WL 

Pi 

Pz 

Pi 

Pi 

0 00 

1 330000 

1 330000 

t 330000 

1 330000 

0 25 

2.657353 

1 856834 

1 856834 

2 657353 

0 50 

2 370394 

1 949367 

1 949367 

2 370394 

0 75 

2 323569 

1 978622 

1 978622 

2 323569 

1 00 

2 309238 

1 990668 

1 990668 

2 309238 

1 25 

2.303136 

1 996640 

1 996640 

2 303 1 36 

1 50 

2 300000 

2 000000 

2 000000 

2 300000 

. 


TABLE  II 

SECTION  LINES  AT  SEVEN  STATIONS 


Station 

a 

P 

7 

<5 

FP 

1 200000 

4.000000 

4.000000 

0.000000 

3} 

1.199377 

1 892324 

1.892324 

0.000000 

65 

1 197513 

0.4801 15 

0. 480115 

0 000000 

10 

1 188297 

0 000000 

0.000000 

0.000000 

13| 

1 042503 

0.191434 

0.191434 

1 182753 

16| 

0 645218 

3 023980 

3.023980 

1.547807 

AP 

1 1 

0 950000 

10.000000 

10.000000 

1.280000 

Station 

£ 

A 

1 

B 

C 

1 1 

FP 

0 100000 

0.006000 

0 000000 

0 013000 

35 

0 276175 
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WATERLINE  OFFSETS 
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Figure  2 Comparisons  biiwnn  Tabular  Data  and  Smooth  Simulations  • tabular  offsets 
from  NSRDC  Report  No  \7\2  waterline  offsets  of  the  smooth  simulations 


Figure  1 Section  Lines  of  Smooth  Simulation  raelie  al  r<  present,.! mn  orthonormal  polynomial  representation 
Stations  of  the  section  lines  are  indicated  in  Figure  1 
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DIHEDRAL  FLOW 


DIHEDRAL  FLOW 


That  the  velocity  is  infinite  for  potential  flow  over  a convex  dihedral  angle  can  be 
understood  from  an  analysis  in  terms  of  complex  variables  Let  ic  be  a complex 
potential  which  is  expressed  in  terms  of  its  real  and  imaginary  parts  by  the  equation 

w = u + iv  ( 1 ) 


Let  z be  a complex  position  which  is  expressed  in  terms  of  its  real  and  imaginary 
parts  by  the  equation 

z = x + iy  (2) 

If  w is  an  analytic  function  of  z.  then  dw/dz  is  independent  of  the  direction  of  dz. 
and  u.v  satisfy  the  Cauchy-Riemann  equations, 

du  dv  dv  du 

--  + I - — - i t-  (3) 

dx  dx  dy  dy 


Differentiation  shows  that  u.  v satisfy  the  Laplace  equations, 

d2u  d2u  d2t  d2i 

dx 2 * dy 2 ~ ° dx2  ' dy2 


( 1) 


If  dz  is  a line  element  in  the  z plane,  then  the  flux  of  Xu  across  a line  with  normal 
i dz  is  the  change  in  v along  the  line  element  dz  Thus  if  u is  a velocity  potential 
such  that  velocity  is  the  negative  gradient  of  u.  then  v is  the  stream  function  such 
that  the  flux  between  streamlines  is  the  difference  in  i 
Let  w be  expressed  by  the  equation 

w = i zy  (5) 


where  7 is  a noninteger  exponent  By  Euler's  theorem  the  complex  parts  are  connected 
by  the  equation 

u * iv  r^sin  y<f>  - 1 r’cos  y<{>  (6) 

where  r.  <f>  are  cylindrical  polar  coordinates  such  that 

1 r cos  0 y r sin  0 (7) 

then  the  stream  function  is  zero  where 

70  (8) 

For  a convex  dihedral  angle,  the  range  of  0 is  more  than  and  7 is  less  than  unity 
Since  derivatives  with  respect  to  r contain  the  factor  ry  the  derivatives  become 
infinite  as  r • 0 

It  is  not  necessary  for  the  fluid  to  adhere  to  the  dihedral  boundary  in  an  ideal 
fluid  In  a more  realistic  flow  there  would  be  separation  with  the  origination  of  a 
vortex  sheet  from  the  vertex  of  the  dihedral  angle  Imagine  what  a network  of  vortex- 
sheets  would  trail  behind  a polyhedron' 
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SPLINE  CURVE 


The  shape  of  the  non  linear  spline  is  computed  with  finite  difference  techniques  by 
Malcolm10.  The  shape  of  the  physical  spline  is  given  preferably  by  elliptic  integrals  for 
which  there  are  efficient  subroutines. 

The  spline  is  a long  slender  elastic  bar.  It  is  bent  into  a curve  by  the  application 
of  forces  and  torques  to  its  ends.  Let  x.  y be  the  Cartesian  coordinates  of  a point  on 
the  curve.  The  shape  of  the  curve  is  determined  by  the  equilibrium  of  forces  and 
torques  on  each  differential  element  of  the  bar.  Let  the  curve  be  divided  into  segments 
of  length  dl  Equilibrium  requires  that  the  net  force  on  any  segment  be  zero,  while 
equality  of  action  and  reaction  requires  that  the  forces  across  a section  shall  be  equal 
in  magnitude  but  opposite  in  sign.  Let  / be  the  force  which  one  side  of  the  bar  exerts 
on  the  segment  of  length  dl.  Then  / also  is  the  force  which  the  segment  of  length  dl 
exerts  on  the  other  side  of  the  bar 

Let  the  constant  force  / be  directed  along  the  y-axis.  The  force  / exerts  a torque 
of  magnitude  f dx  on  the  segment  of  length  dl.  The  torque  must  be  balanced  by  a 
difference  dM  in  the  moments  which  are  applied  to  the  ends  of  the  segment  The 
moment  M is  proportional  to  the  curvature  of  the  bar. 

If  no  force  is  applied  to  the  bar,  then  the  moment  M is  constant  The  bar  is  bent 
into  an  arc  of  a circle. 

If  a force  is  applied  to  the  bar,  then  the  moment  U is  a linear  function  of  x Let 
the  origin  of  x be  located  at  a point  where  M is  zero  Then  the  moment  M is  given  by 
the  equation 


M - fx 


The  moment  U is  given  in  terms  of  the  curvature  of  the  bar  by  the  equation 


M 


K 


de 

dl 


(1) 

(3) 


where  K is  the  bending  modulus  of  the  bar,  and  0 is  the  angle  which  the  tangent  to 
the  curve  makes  with  the  x axis 

The  angle  0 is  defined  by  the  equation 


0 


tan 


riy 

dx 


Differentiation  with  respect  to  x leads  to  the  equations 

dzy 


dO  dxz 


(3) 


(4) 


(5) 


1 


while  the  radius  of  curvature  R is  given  by  the  well-known  equation 

dzy 

1 dd  dx2 

R dl 


(i)  r 


The  spline  curve  is  determined  therefore  by  the  equation 

dy 

dxz  f 


( \ dx  / ) 


, — x 

l K 


Integration  leads  to  the  equation 


dy 

dx 


\ 


m 


where  a is  an  arbitrary  constant  and  c is  defined  by  the  equation 

/ 

c = — 

2A' 

Solution  for  dy  dx  and  integration  lead  to  the  equation 

* f . 

J \ 1 a cx2  \ 1 - a + cx2 


The  substitutions 


and 


(V'f 


COS  0 


k 


(V) 


y F(<t>.  k)  ^'E(tt>.k) 


(6) 


(7) 


(8) 


(9) 


(10) 


(ID 


(12) 


onvert  the  integration  into  elliptic  integrals  as  expressed  by  the  equation 

1 

v 2c 

where  the  elliptic  integrals  of  the  first  and  second  kinds  are  defined  by  the  equations 


(13) 


F(<(>.k) 

f*  d * 

Jo  x 1 t-2sin “v 

A'(0. k) 

\ 1 fc2sin2^  dy: 

Jo 

(14) 

(15) 


The  elliptic  integrals  give  the  locus  of  a free  elastic  spline  which  is  subject  to  forces 
alone  at  its  ends 
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CALCULUS  OF  VARIATIONS 


Let  x be  an  independent  variable  and  let  y,  y' , y"  be  a dependent  variable  and  its 
derivatives  with  respect  to  x Let  f(x,y,  y' , y")  be  the  integrand  of  the  integral 

f fix.y.y’.y")  dx  (16) 

J a 

where  a,  b are  fixed  limits  of  integration  Let  the  integral  be  stationary  with  respect 
to  any  arbitrary  variation  as  expressed  by  the  equation 

<5  I f(x,  y,  y',  y")  dx  = 0 (17) 


Variation  in  the  integral  is  derived  from  variation  in  y.  y . y"  as  expressed  by  the 
substitutions 

y - y + 6y  (18) 


y'  - y'  + -j  <5y 

dx 


y"  - y"  + , 2 


Substitution  in  the  integral  leads  to  the  equation 

<5  | f(x,y.y',y")  dx  | + - 


df  df  , ) 

°y  + oy  + TT.  °V  \ a* 

ay  dy  ) 


and  integration  by  parts  leads  to  the  equation 

A*.  ».*■.»")*  [ {(g)-  i‘(iy)\Syd‘ 

■o:  £(%)«%)«■:  «* 

Thus  if  dy,  6y'  are  zero  at  a.  6 then  the  integral  is  stationary  for  any  arbitrary  <5y  only 
if  / is  a solution  of  the  Euler  differential  equation 

(df)  *(.V)  0 (23) 

\ dy  I dx \dy' I dx2  \ dy  " } 

Solution  of  the  differential  equation  determines  the  relationship  between  x and  y for 
a stationary  integral 


MINIMUM  CURVATURE 

In  the  case  of  the  integration  of  the  square  of  the  curvature  along  the  length  of 
a curve  the  integral  is 

fb  (y")2 

V & ’ (Ot) 


r ax 

ll  + (y)2!2 


3 


1 


Inasmuch  as  the  integrand  does  not  depend  upon  y,  the  Euler  equation  can  be  integrated 
to  give  the  equation 


a/  d 

— + — , 

dy'  dx  \ dy 


Substitution  leads  to  the  equation 

2 y'" 


V-  = 

dy"  I 

5y'(y")z 


constant 


1 1 + (y')2i2  f i • (y')2l2 


constant 


(25) 


(26) 


If  this  equation  is  multiplied  throughout  by  y",  it  may  be  integrated  to  give  t In- 
equation 


(y")2 

U + (y')2i3 


A ■ By 

\ i • (Vr 


(27) 


where  A.B  are  arbitrary  constants  Let  the  coordinate  axes  be  oriented  so  that  th< 
x axis  is  in  the  direction  of  the  vector  .-li  - By  Let  the  angle  0 be  the  angle  between 
the  tangent  of  the  curve  and  the  x axis  Then  the  angle  0 satisfies  the  equation 


/ ^0  \2 

C sin  0 

' dl  I 


(28) 


where  C is  the  magnitude  of  the  vector  Ai  + By  The  x coordinate  is  related  to  the 
distance  l along  the  curve  by  the  equation 


dx 

dl 


cos  6 


Substitution  and  integration  leads  to  the  equation 

i i 

2 sin20  = Czx 


(29) 


(80) 


where  the  origin  of  x is  selected  to  be  that  point  where  0 is  zero  Squaring  of  both 
sides  of  the  equation  leads  to  the  equation 


sin  0 


y 

\ 1 ■ (y ')2 


(81) 


Comparison  of  this  equation  with  the  equations  for  the  physical  spline  shows  that  t In- 
curve of  minimum  square  of  the  curvature  is  a special  case  of  the  elastic  curve 
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